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We study the consequences of CPT and lepton number violation in neutrino sector. For CPT 
violation we take gravity with which neutrino and antineutrino couple differently. Gravity mixes 
neutrino and antineutrino in an unequal ratio to give two mass eigenstates. Lepton number violation 
interaction together with CPT violation gives rise to neutrino-antineutrino oscillation. Subsequently, 
we study the neutrino flavor mixing and oscillation under the influence of gravity. It is found that 
gravity changes flavor oscillation significantly which influences the relative abundance of different 
flavors in present universe. We show that the neutrinoless double beta decay rate is modified due 
to presence of gravity— the origin of CPT violation, as the mass of the flavor state is modified. 
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I. INTRODUCTION 



The oscillations between different kinds of neutrino and antineutrino flavor have been observed form solar, atmo- 
spheric and LSND data. The three pieces of observation indicate three values of mass-squared difference of three 
different orders. With three families of neutrino, one can obtain only two independent mass-squared differences. 
Therefore, observations require the introduction of fourth neutrino which must be sterile in the standard model. But 
many difficulties arise with the introduction of fourth neutrino as discussed in literature (e.g. see 

As an alternate proposal to accommodate the results, many authors have proposed CPT violation in neutrino 
sector 0, H, 0] • One can either introduce a new particle (sterile neutrino) or allow CPT violation to take care of all 
experimental results with present data. However, very recently, MiniBooNE results have been declared which shows 
that LSND experimental results can not be explained simply by neutrino oscillation. Hence, it raises many other 
questions to be answered @. 

Previously, many authors have explored the consequences of CPT violation in neutrino sector (e.g. see @, 0|)- The 
nature of mass of neutrino has been studied [?} if CPT is violated. It has been shown that if CPT is violated then 
neutrinos no longer remain Majorana particles even if they have a Majorana mass i.e. they violate lepton number. 
But the actual physical condition which can lead to CPT violation was not discussed satisfactorily. Proper situations 
in which CPT can be violated were discussed later by many authors [1, H, [l(| HH, EE 03 • They have pointed out 
that CPT violation in neutrino sector can occur due to the spin-gravity coupling. 

In the present paper, we plan to obtain the mass matrix for neutrino sector with Majorana mass in the presence of 
gravity. It has already been seen that spin gravity coupling has an extra contribution to the effective mass of neutrino 
and antineutrino. The Majorana type mass of neutrino alone can cause lepton number violation i.e. the mixing of 
neutrino and antineutrino states. Without an interaction with gravity, the neutrino-antineutrino mixing angle is 7r/4. 
However, in presence of background gravitational field, mixing angle changes and depends on the strength of CPT 
violation. In a similar fashion, flavor oscillation is also altered by gravitational effect. As a consequence of mixing, 
the neutrino-antineutrino oscillation can take place only in presence of gravity. We also show that gravity affects the 
neutrinoless double beta decay rate even if two neutrino flavor mixing is taken. 

We organize the paper in the following manner. In the next section, we review the CPT violating nature of the spin- 
gravity interaction for neutrino. In the section Hill we obtain the mass matrix for neutrino traveling in a background 
gravitational field. From that we find the neutrino-antineutrino mixing and oscillation. Subsequently, we investigate 
effect of gravity on flavor oscillation and its application to that in early universe and around black holes in the section 
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IIV1 Applying this mixing, we recalculate the neutrinoless double beta decay rate in the next section. Finally, we 
summarize what we obtain and discuss their implications in the section fVTl 



II. COUPLING TO CURVATURE 



The CPT violation mechanism due to the spin-gravity coupling of fermions has been discussed earlier [8L]9, 1C , 13 1. 
How differently gravity interacts with neutrino than antineutrino has been shown by a detailed calculation [10 . Il3l . 1 14 1 . 
For completeness, here we revisit the idea very briefly. 

The general invariant coupling of spin-1/2 particles to gravity is described by the Lagrangian 

S3, S3 



where the covariant derivative is 



and the spin-connections are 



Here 



£ = x -<j ( ^* 7 a Da * - *m*J , 



D a = [ a — -U>bcaCT bC 



Ubc 



ebx (d a e 



(1) 



(2) 



(3) 



(4) 



All the above equations have been written in a local inertial frame which is flat along the entire geodesic. The Roman 
alphabets denote the flat space indices, while the Greek alphabets denote the curved space indices. Here e-s are 
the vierbeins connecting curved and locally flat spaces and obey the relations e^e l/a - 
represents the inertial frame Minkowski metric and gV the curved spacetime metric. 
Thus the Lagrangian can be rewritten as 



a afi e b 



T) , where r\ 



a b 



C = det(e) 4- ( - 7 a da -m + 7°7 # 



with 



T>d ^abcd, . 



(5) 



(6) 



In a local inertial frame, the effect of gravity on fermion is only an axial vector interaction term involving the 
gravitational field B a given by Eq. ©. The Lagrangian given by Eq. ([5]) has two parts— the free part and the 
interaction part. The interaction part is a coupling to the field B a which is constant in a local inertial frame because 
that arises from the background gravitational field. 

If B a is constant in a local inertial frame, then the interaction term violates CPT as well as particle Lorentz 
symmetry, although it is invariant under observer Lorentz transformation (see e.g. [HI)- It has been shown that only 
a special form of the background gravitational field can give rise to such CPT violating interaction term. For example, 
around the rotating black holes [9. Hcj|. or in anisotropic early universe 13], a suitable background field exists which 
causes CPT violation in interaction with spin-1/2 fermions. 

The interaction term in the Lagrangian given by Eq. ([5]) involving B a contains 7 . If we consider the standard 
model neutrino with the Majorana mass term, then right-handed neutrino and left handed antineutrino are absent. 
This leads the gravitational interaction to — ^l^^l and ^ c L ^ a '^ c L for neutrino and antineutrino respectively. Here 
c superscripted ^ represents the charge conjugated spinor or the spinor for antiparticle and the subscript L denotes 
the left handed component. Consequently, the dispersion relations for neutrino and antineutrino become different 
respectively given as 



E » = \/(p-B) 2 +m 2 + B , 
E v . = J {p + B) 2 + m 2 - B 



(7a) 
(7b) 



with momentum p. 
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III. MASS MATRIX 

If we consider Majorana neutrino, then in Weyl representation neutrino spinor can be written as 

•-(£)■ (8) 

where ^>£ and tpL are two component spinors for antineutrino and neutrino respectively which are lepton number 
eigenstates with eigenvalues -1 and +1 respectively. Here we assume that neutrino is a left handed particle only. 
Hence, hereafter we shall omit the subscript L. 

In terms of two component spinors, the Majorana mass term can be written as 

*<m*= (^^)(_°„- m )(t; c ) 

= — Tp^mip — tfj^mip . (9) 
Now in gravitational field the Lagrangian density can be written as 

(- 9 )-" 2 £ = (r< *«)f>V5,,(£) + («" «t)T'Bc(*° 

- (*- *')rVS,(!f) + (^+*)(*_° ab )(t: 

- (*-**•) %) w 

where 

= (a ,ft+ 7 5 B l ). (11) 

Therefore, in Lagrangian density we obtain terms containing Bo4> cj< ip c and Botp^tp. These two terms do not vio- 
late lepton number. Hence these terms can be looked as lepton number non-violating mass of Majorana neutrino. 
Therefore, writing the mass terms together we obtain the Lagrangian density 

(-,)-«»£ = ( ) • yy ** *> ) ( ~Bq ){%) < 12 » 

Hence, in background gravitational field the Euler-Lagrange equation for neutrino and antineutrino is 

- (^x)O o - <i3 » 

Therefore, the mass matrix of the neutrino-antineutrino sector is given by 

(14) 

Hence we see that if one assumes neutrino to have solely the Majorana type masses, then it acquires lepton number 
non-violating type masses, equal but opposite in sign, while propagating in a gravitational field. Moreover, we see 
that in this case ij> no longer remains a mass eigenstate. 

A. Neutrino-antineutrino mixing and oscillation 



The mass matrix M., given by Eq. (|14p , is Hermitian and can be diagonalized by unitary transformation. Then 
the mass eigenstates will be some admixture of ip and tp c . Hence, neutrino and antineutrino states couple together. 
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We find these two mass eigenstates v\ and v 2 as 



X - | (b + y/F* + m^j m + m|V)| (15a) 



W v 

M = ^\ v) + [B„ v '.^- :i " :J ) '•}}. 'ir,h = 



So + JbJ+^A \i>) 



with eigenvalues 



Here 



mi,2 =Ti/B 2 +m 2 . (16) 



iV = \j2Bl + 2m 2 + 2B yJ Bl+m 2 (17) 

is the normalization factor. In a more convenient way, one can write 

= cos(9 \i) c ) + sm9 \ip) (18a) 
\v 2 ) = -sm9\iP c ) + cosd |V) (18b) 

with 

tanfl = , (19) 

Bo + y/B* + ™ 2 

Then one can write \ip c ) and as the superposition states of and \v 2 ) in the following manner 

\tp c ) = cos 9 - sin6» |z/ 2 ) (20a) 
= sin0 + cos6» |i/ 2 ). (20b) 

Basically, v\ — — ia 2 v 2 as ?/> c = —i<J 2 i\). Therefore, !/i and 1/2 together actually describe a single Majorana particle 
as ip and V c do. Then, if we construct a new four component neutrino spinor as 



(21) 



it retains Majorana nature, i.e. v c = v. This is obvious as the new spinor v is only a transformed spinor from the 
previous one '5 and the old and new spinors are related by a unitary transformation 



cos V sin ( 
— sin 9 cos ( 



(22) 



Although, v\ and v 2 are not lepton number eigenstates, as evident from Eq. (|15|) . the mass terms in the corresponding 
Lagrangian of V\ and v 2 are lepton number conserving. In absence of gravitational field, neutrino and antincutrino 
mix in the same angle. 

Hence, if there is any lepton number violating interaction, then we obtain two neutrino mass eigenstates v\ and 
v 2 which are superpositions of neutrino and antineutrino states given by Eq. (|18|) . As energies of neutrino and 
antineutrino in gravitational field are different, there will be an oscillation between v\ and v 2 . At any time t, the 
oscillation probability is given by 

V(t) = sin 2 29 sin 2 <5(i) (23) 

where 

S(t) = l^-^l* . (24) 
when E v and E v c are given by Eq. ([7]). Hence, for ultra relativistic neutrinos by putting the value of we obtain 

V ® = 77T^ sm 2 {(B -\B\)t} (25) 
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for neutrino and antineutrino of same rest mass to. Thus the oscillation length is 

A= (26) 

Bo-\B\ 

This depends only on the strength of the gravitational field. If we consider neutrinos to be coming out off the inner 
accretion disk around a spinning black hole of mass M = 10M©, then Bq — \B\ = B = 10~ 19 GeV jl4j| . which leads 
to A ~ 10km. If the disk is around a supermassive black hole of M = 1O 8 M in an AGN, then A may increase to 
10 10 km, depending upon the size of inner edge where from neutrinos come out and angular momentum of the black 
hole. Therefore, an oscillation may complete from a few factor to hundred Schwarzschild radii in the disk producing 
copious antineutrino over neutrino and may cause overabundance of neutron and positron. However, neutrinos around 
a primordial black hole of mass M p [38| could lead to an oscillation length as small as A ~ 10~ 16 km < 100M P for 
r < lOOMp. 

Therefore, from Eqs. (p~8|) and (| 1Q[) we see that for nonzero value of Majorana mass (m ^ 0), neutrino and 
antineutrino combine to give two new states V\ and v%, mass eigenstatcs with different mass eigenvalues. Therefore, 
in these two states neutrino and antineutrino are mixed. In the absence of gravity, the mixing angle is 7r/4 which is 
evident from the Eq. (|19[) . On the other hand, if m = 0, i.e. there is no lepton number violating interaction, then 
neutrino and antineutrino do not couple at all [see Eq. (|13[) ]. This means two component neutrino and antineutrino 
described in Eq. © themselves are mass eigenstates. In this case, presence of gravity which is CPT violating, splits 
these two eigenstates with two different mass eigenvalues. 

We see that although initially gravitational interaction and the Majorana mass term explicitly have different effects, 
one to violate CPT and another to violate lepton number, both of them contribute in the same manner to split the 
mass eigenstates (see Eq. (TTBIO . Moreover, in gravitational field, since neutrino and antineutrino acquire different 
effective masses, gravitational field coupled to neutrino spin may have some lepton number violating nature implicitly. 
This has been illustrated in literature [l(| [Hj]. We also see from Eq. (|25|) that presence of gravity leads to oscillation. 
Without gravitational field, the lepton number violating interaction alone can not cause this oscillation. 

B. Oscillation with lepton number conserving mass 

For Majorana neutrino when the sterile components are neglected, the Dirac mass is of no relevance. This is exactly 
the case we are considering in our present paper. We have also seen that the gravity can induce an effective mass 
for Majorana neutrino which does not violate lepton number. If we consider the diagonal term of the mass matrix of 
Eq. ^ to be non-zero, then we obtain an extra mass term which conserves the lepton number. With this mass term 
included, the effective mass matrix will take the form 

where m n is the so called lepton number non- violating mass. 

With this mass matrix, the mixing angle of neutrino and antineutrino does not alter. However, the mass eigenvalues 
of the mass eigenstates become 

m n (i,2) = m n T \J Bl +m 2 . (28) 

In this case, we can consider oscillation between neutrino and antineutrino as they are combination of v\ and V2 
which evolve differently with time. For this oscillation, the oscillation probability at any time t is given by 

T n {t) = sin 2 26> sin 2 e(t) (29) 

where 

<t) = l^M, (30) 

when E\ and E2 are energies of two neutrino mass eigenstates with momentum p. In the ultra-relativistic limit 
(\p\ » m) 

2 

£ ^ = i*i +! ^gr (31) 
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Again assuming |p| ~ E 



m n2 - m nl _ 2m„ y/B 2 + 171 2 



E2-Bk= nZ 2E nL = ■ (32) 



Hence the the oscillation probability at any time t is given by 



and the oscillation length is given by 



(34) 



Hence we see that this neutrino-antineutrino oscillation probability does not depend on the spatial part of the 
gravitational field (gravitational vector potential) B but depends on the temporal part (gravitational scalar potential) 
-Bo only. However, there is nothing new in it. It was already shown that neutrino asymmetry and then leptogenesis in 
early universe arises due to non-zero Bq [131 ] independent of B. The gravitational scalar potential Bq is non-zero if the 
background spacetime has at least one nonzero off-diagonal spatial component when the set of coordinate variables is 
{t, x, y, z}. Existence of such a component in a spacetime may be due to anisotropy which is the case for the Bianchi 
model. Moreover, off-diagonal spatial components also govern in presence of primordial quantum fluctuations in the 
Robertson- Walker spacetime. This is basically the tensor perturbation to early universe. Therefore, it is confirmed 
that the gravity induced leptogenesis is possible only if the gravitational scalar potential is non-zero. 



IV. EFFECT OF GRAVITY ON FLAVOR MIXING AND OSCILLATION 

A. Flavor mixing 

In previous section we have seen that interaction of neutrino and antineutrino with gravity gives rise to neutrino 
mass eigenstates which are superposition of neutrino and antineutrino states given by 



= UlJfw (35) 



where 



14 „ = ( C0S ^ - Sin /^V (36) 
e ^ ^ sin e , M cos 6 e ^ J K ' 

It is clear from our above discussions that z/ e „ here are four component spinors. In other way, it can be stated that 
in gravitational field neutrino state \P of mass m is modified to the state v of mass yj Bq + m 2 as described in section 
IIIII Therefore, recalling Eq. (|2Tj) we construct flavor eigenstates v e and under gravity for electron and muon 
neutrino respectively whose components are v e \, u e 2, v^i and with masses m e i, m e 2, m M i and m M 2 respectively. 
In terms of gravitational coupling these mass eigenvalues are expressed as 



m ( e , M )2 = yB$+m 2 e ^ (37) 

where m e , are the Majorana masses for electron and muon neutrino respectively, analogous to m of section III. 
The corresponding mixing parameters are given by 

tan0e, M = me ' M (38) 

B + ^B 2 + m 2 ^ 

Now we consider the corresponding two flavor mixing. The states v e and are coupled by a Majorana mass term 
m e/J . Then the mass term in the Lagrangian density is considered as 
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- v^m^Vei - vl-ym^v^x + v^ 2 m e ^v e2 + v^m^v^ . (39) 

Here we assume that the Majorana mass matrix, coupling v e and v^, is Hermitian and diagonal. We also assume that 
m eM in v^v^i and v\ 2 v^ 2 are same for computational simplicity, while our main goal is to investigate any gravity effect. 
As our main aim is to study the effect of curvature to the oscillation phase, even for the convenience of transparent 
understanding of sole effect of gravity we prefer to consider the nongravitating part as simple as possible. This mass 
term gives rise to two mass matrices M.f (1,2) given by 

M = ( m el -m efl \ M ( ^ ™w ) (40 ) 

which mix v e \ {v e2 ) and (1^2) • The mixing, leads to two mass eigenstates f\\ (/12) and /21 (/22), is expressed as 

/ = ftvf (41) 



where 



and T is given by 



= cos*,., . (43) 

\_ Sin 0i,2 COS 01,2 / 

Here the subscript 1 and 2 refer to the flavor mixing between mass eigenstates subscripted by 1 and 2 respectively. 
Hence 

tan 0i, 2 = T2 (44) 

(«>e(l,2) - to m(1,2)) + J K(1,J) - «V(i,2)) 2 + 4m^ 

Thus we obtain all together four mass eigenstates xii X27 X3 an d X4 described as 

with mass eigenvalues Mi, M 2 , M 3 and M4 given by 

Mi, 2 = m ( i, 2 )i = - |(m e i + m M i) ± y/ (m e i - rn^i) 2 + 4m^| 

M 3 , 4 = m(i, 2 )2 = X {( m e2 + "V*) ± y/ (m e2 - m^) 2 + 4m^| . (46) 

Now from Eq. ([37]) 



m el +m^ = - [jB 2 + m 2 e + jB 2 +m 2 ) = — (Y e + Y^) 



m e2 + m^ l2 = [ \ B$ + m 2 e + J B$ + rn 2 ) = Y e + Y f 



m e i - nfyi, = - \/ B^+m 2 e + JB^ + m 2 = -Y e + 



m e2 - m M 2 = J B 2 + m 2 e - J B 2 + m 2 = Y e -Y tl (47) 
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where 



If we define 



then 



Then we have 



Y=JB 2 + m 2 . (48) 



2m a = + Y e , 2m t = — Y e (49) 



m e i + rrifii = -2m a , m e i - m M i = 2m ll 

m e2 + m M 2 = 2m a , m e2 - m M2 = -2mj. (50) 



Mi = - m a + \/mf+m^, M 2 = - m a - Jm^+m 2 ^, 



M 3 = m + Jmf + m 2 M 4 = m Q - Jmf + m 2 (51) 



and 



taa^i ia = (52) 

Note that at -Bo = all the above results reduce to that of conventional two flavor mixing without gravity effect. For 
example, at i?o = 

V'e./i = -^=(l>(e,At)l +' / (e,/i)2)) 

^e, M = -^0(e, M )l - ^(e,M)2)- (53) 

Then substituting f-s by -0-s in Eq. (|39[) with TO( e , M )i = — ^e,^ and m^ e ^ 2 = ^e,/i one obtains the standard mass 
Lagrangian density with spinors ^ and V' [25l | 

£ m = --(^'fflei + ipfme^ipu, + %p L ^m efl ip e + i/^m^W) + ^- c - ( 54 ) 

The readers compa ring this result with that in 25] should not confuse with difference in notation in the present 
paper from that in [25j. Here v e ^ are the Majorana neutrino fields of electron and muon typed under gravity whose 
oscillation would be interesting to study. 

The most general flavor mixing mass matrix in presence of gravity is given by 

-B I -M 

Jo: 

when I is the 2x2 unit matrix and 



M ^=\ — M B„I ) ' ^) 



M = I m, I = Ue ' diag(mi < m2) ' U «' (56) 



where TJg and m\ 2 are the mixing matrix in vacuum and the neutrino masses respectively in absence of curvature. 
The the corresponding flavor state is given by 



I % 



■ (57) 



However, for the sake of simplicity we take an assumption such that the 4x4 mass matrix is resolved into 2x2 
block form. The motivation is to study the gravity effect in a transparent manner. In order to do that we consider 



9 



the mixing part by part as described above: first we take the gravity effect on the flavor eigenstates and then we mix 
them up. The argument for this assumption is as follows. We detect neutrinos by weak interactions. In standard 
experiments, we detect and produce them via charge-current interaction. As the flavor eigenstates take part in 
charge-current interaction, we detect and produce them only. In our calculation, the original flavor eigenstates are ip 
and Tp c or as a whole ^. Therefore, we plan to discuss the oscillation probability and corresponding length among 
states vp. However, when we consider them under strong gravity, they couple and no longer have been a definite 
mass state, rather modify to a state v. Hence it is worthwhile to study the oscillation properties for states v as well. 
Indeed the solution of the Dirac equation in presence of gravitational interaction is different than that in absence of 
gravity. Hence while performing the oscillation experiment under gravity, it can be assumed that gravity will affect 
the initially produced flavor states and in between the detection and production the flavor mixing Hamiltonian will 
act upon these gravitationally modified flavor states. Therefore, even though they are produced as flavor eigenstates 
\l/ via weak interaction, if system is under strong gravity, then they are modified to v, as described in §111. A, B, which 
may be observable in presence of strong curvature. Indeed v = "J when Bq = 0. If we put Bq = in our oscillation 
probability and oscillation length, described in §IV.B below, then we obtain the standard expressions of those in flat 
space. This validates our assumption of mixing scheme and verifies that our results for v are physically same as that 
for ^ under gravity, as will be discussed in §IV.C. 



B. Flavor oscillation 



The sets of two Majorana neutrino flavor eigenstates are described as 

\v e \) = cos 0i | xi) - sin0i|x 2 ) 

|f M i) = sin0i|xi) + cos 0i |x 2 ) 

and 

\Ve2) = COS 02 |X3) ~ Sm0 2 |x4) 

| I/ M 2) = Sin 02 |X3) + COS 02 |X4)- 

Let us now consider the states given by Eqs. ([55]) . The oscillation probability of those states is 

Vfgi = sin 2 20i sin 2 6f g i (t) , 

where 



Sf 9 i(t) 



\M\-Ml\ 



t 



and E is the energy of the mass eigenstates. Hence the oscillation length is 

AirE 



\M\-Ml\ 



Similarly, for the states of Eqs. (|59"|) the oscillation probability is 

V fg2 = sin 2 202sin 2 5 /s2 (X): 

where 



and the oscillation length is 



_ \M!-Mj\ 
SfgiW ~ 1£ * 

AttE 



Vfl2 



\MI-Ml\ 



Now from Eqs. (H7J), (J49J) and (JHTj) we obtain 



AM 2 = \M\ - M|| = \M\-M\\ 



Bl + m 2 



Bl + ml 



1 



Bl + ml 



Bl + ml 



(58a) 
(58b) 

(59a) 
(59b) 

(60) 
(61) 

(62) 
(63) 
(64) 

(65) 



(66) 
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From Eq. (|52[) , we obtain 

2 

sin 2 20i = sin 2 20 2 = sin 2 20 = ^ _ (67) 

so that the oscillation probability and the oscillation length are same for the two cases respectively given by 

V fg = sin 2 20 sin 2 6 fg (t), (68) 

Sg = AAV' (69) 

However, Vt-s are related to f-s. Hence once we obtain the oscillation probability and length between v e i (or v e2 ) 
and (or v^), the results can be converted to that for ip e and Vv given in the next subsection. 



C. Oscillation probability between initial flavors produced via weak interaction 

It is evident from Eqs. (f3"5| and (j4"Tj) that we can express the initial flavor states f e and ^ in terms of mass 
eigenstates x~ s given by 

= COs9 e COS 01 Xl ~ COs0 e Sm 01 X2 — Sin0 e COS 02 X3 + Sin# e Sin 02 X4j 

0£ = cos0 M sin 0i xi + cos0 M cos0i X2 - sin0 M sin0 2 X3 ~ sin0 M cos0 2 X4, 
ip e = sm.6 e cos 0i xi — sin 9 e sin0i X2 + cos(9 e cos 02 X3 ~ cos(9 e sin 02 X4j 
0^ = sin# p sin 0i xi + sin0 M cos0i X2 + cosd^ sin 02 X3 + cos# p cos0 2 X4< (70) 

In short, we can write this as 







(xi 






X2 


X3 






\ Xi 



(71) 



where 



/ cos 9 e cos 0i — cos 9 e sin tj>\ — sin 9 e cos 02 sin 9 e sin 2 



cosf 



sm (pi cos On cos (pi 



sm On sm ip2 



sm fe cos 01 



sm e sm 0i cos 9 e cos 02 



T = 

\ sm sm 0i sm 0^ cos 0i cos f M sm 2 
From this the particle part ^ can be written as 

i'e = Tel Xl + T e2 X2 + T e3 X3 + T e4 X4, 

^fl = T ttl Xl + T li2 X2 + 7^3 X3 + 7^4 X4- 

Hence we compute the oscillation probability between tp e and tp^ 

4 



— sm cos <p2 

— COS # e Sin 02 
COS Ofj, COS 02 



where 



Sij(t) = 



i<j=l 



I AM 2 



AE 



AM 2 = Mi 



Mi 



(72) 



This resembles the oscillation probability one obtains in the case of 3 flavors mixing or in general the case of ./V 
mixing. With actual expressions of T^^j-s from Eq. (|70p we obtain 

Vi g = sin 2 8 e sin 2 9^ sin 2 20i sin 2 5i2(i) + cos 2 9 e cos 2 9^ sin 2 202 sva 2 5^i(t) 

— isin2# e sin 20^ sin20i sin20 2 {sin 2 5ia(t) — sin 2 5u{t) — sin 2 5 2 3* + sin 2 Suit)} . 



(73) 
(74) 

(75) 
flavors 

(76) 
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From Eq. ([51 



and 



\M 2 -M 2 \ = \M\-Ml\ = AM 2 = Am a Jm 2 +m^ (77) 



\M 2 -M 2 \ = |Mf-Mf| = 0. (78) 

Moreover, from Eq. (|52|) 

sin 20i = —sin 202 = sin 20 (say). (79) 
Hence the oscillation probability is given by 



V ig = sin 2 20 <^ (sin 6 e sin M + cos e cos 6» M ) 2 sin 2 ( -f^-i ) [• . (80) 



/ AM 



Therefore, we hnd that the oscillation length between ip states is same as that in the case of oscillation between V\ 
(and V2) states: equivalent to each other. 

Now the oscillation length would change for a distant observer due to gravitational redshift. Therefore, the redshifted 
oscillation length 

A/, = ^L- (81) 
V9tt 



If we consider the neutrino flavor oscillation in a rotating black hole spacetime of mass M and specific angular 
momentum a, then at a point (r, 8) 

P 

with 

p 2 = r 2 +a 2 cos9. (83) 

Therefore, it is evident from the above discussions that under gravity the probability of flavor oscillation changes 
significantly depending on the gravitational strength. 

We now consider Majorana mass of electron neutrino m e ~ 5 x 10~ 3 eV, muon neutrino ~ 6.5 x 10 -3 eV and 
mixing m efi ~ 3.5 x 10~ 3 eV. These values are consistent with the solar neutrino oscillation data. Thus, without 
gravity the oscillation probability between electron and muon neutrino with AM 2 (Bq = 0) = 8.2 x 10~ 5 cV 2 , as from 
solar neutrino data [261 ]. is given by 

, n 9 /8.2 x 1CT 5 eV 2 \ 
?f = VfgiB = 0) ~ 0.956 sin 2 f — t\ . (84) 



The corresponding oscillation length 



4ttE 

A / = A /g (j^O) ^ 2xio _ 5eV2 . (85) 



D. Oscillation around black holes 



We first recall the gravitational scalar potential computed earlier in the Kerr geometry [lj 



B° = -^t (86) 

for a black hole of mass M and specific angular momentum a, where p 2 — 2r 2 + a 2 — x 2 — y 2 — z 2 . If we consider 
neutrinos at around 20 Schwarzschild radius in the spacetime of a primordial black hole of mass 10 22 gm, then from 
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Eq. (|86|) the gravitational field is computed to be ~ 10 2 eV which is comparable to neutrino masses. If we specify 
the gravitational field B ~ 5 x 10~ 2 eV, then from Eqs. ([55)) the oscillation probability becomes 

V f BH ^ 0.999 sin 2 f ? * ^ ^ t) . (87) 



The oscillation length, from Eq. (|69|) . is determined as 



A/BH " TxlO^eV^ (88) 

Hence, the oscillation length decreases by almost an order of magnitude. 

If Bq is much higher than neutrino Majorana masses, then from the Eq. 1|67|) we have sin 2 2<j) ~ 1 and AM 2 ~ 
4m eAt _Bo. Therefore, the oscillation probability 

. 2 /1.4 x 10- 2 cV5 \ 
V f BH = sin 2 ( — °- t ) (89) 



and the oscillation length 



AttE 

XfBH ~ 1.4 x 10-2 eVB - (90) 

E. Oscillation in early universe 

1. Anisotropic universe 

We recall the gravitational scalar potential in the anisotropic phase of early universe [III, Q3] 

_ 4fl 3 S + 3y 2 fl S 3 -2yS 4 

8i? 4 + 2y 2 i? 2 S 2 ' 1 J 

If we consider radiation dominated era with R(t) — (t/to) 1 / 2 and S(t) as an arbitrary constant > 1, then the above 
potential reduces to 



y \t 

when y is the position coordinate, can not be greater than the size of universe of corresponding era, and to is the 
present age of universe ~ 10 17 sec. At the neutrino decoupling age of universe when t ~ 1 sec, y < 10 20 cm, Bq > 10 -8 
eV. Therefore, if Bq ~ 5 x 10 -2 cV, then from Eqs. (|68|) and ([87]) the oscillation probability is given by 



„„„„ . 2 /7 x 10" 4 eV 2 sec\ 
P/EC/, < 0.999 sin 2 f — J . (93) 

On the other hand, at GUT scale when t ~ 10 -35 sec, y < 10 2 cm, Bq > 10 45 eV >> m e ,m^,m efJi . Therefore, the 
oscillation probability with the minimum possible Bq 

2 /1.4 x 10 8 eV 2 sec\ , x 

VfEUauT < sm 2 f — J (94) 

is entirely controlled by gravitational field. Either of Eqs. (|93| and (|94|) clearly proves that flavor oscillation is severely 
altered by the gravity. 

2. Inflationary era of universe with primordial fluctuations 

We know that during inflation primordial quantum fluctuations of the spacetime is classical and the flat Robertson- 
Walker expanding universe may take the form as [13, HI] 
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ds 2 = (1 + 24>)dt 2 - a{tf 



■^dx'dt + ((1 + 2^)% + hij^j dx l dx j 



(95) 



where and ip are scalar, u) i are vector and hy are the tensor fluctuations of the metric. Of the ten degrees of freedom 
in the metric perturbations only six are independent and the remaining four can be set to zero by suitable gauge 
choice. For our application we must have atleast one nonzero gij, when i ^ j = 1, 2, 3, of the metric [To| and thus 
need only the tensor perturbations and we choose the transverse-traceless (TT) gauge h\ = 0, d l hij = for the tensor 
perturbations. In the TT gauge the above perturbed Robertson- Walker metric can be expressed as 



ds 2 



(l + 20) 



dt 2 - a(t) 2 



i(t) 



dx l dt + (1 + 2-0 - h + )dx\ + (1 + 2-0 + h+)dx\ + 



2h x dx\dx2 + (1 + 2ip)dx\ 



(96) 



Therefore, following general expression [T(| [HJ given for any metric the gravitational scalar potential computed for 
this spacetime is 



b° = d 3 h x =d z h x . 



(97) 



The gravitational scalar potential Bq can be expressed as a fluctuation amplitude A x times a wavenumber which 
represents the length scale over which the metric fluctuations vary. The Compton wavelength of the particles in 
the GUT era is much smaller than the average wavelength of the gravitational waves whose wavenumber k ~ H = 
lMgl /2 (T 2 /M Pl ) [13, HI] . H is the Hubble constant at the time of hor izon crossing of mode k. Thus gravitational 
wave background can be considered as a constant amplitude field for the GUT processes. Hence the mean value of 
Bo, as a function of temperature and the primordial tensor wave amplitude A x , can be expressed as 



(B ) =B ~A x k~ A x 1.66 g 



1/2 
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PI 



(98) 



Here = 106.7, is the number of relativistic degrees of freedom, for the standard model. Primordial tensor and 
scalar perturbations contribute to the anisotropy of cosmic microwave background at large angles. The COBE DMR 
measurement [29] of temperature anisotropy AT = SOfiK sets an upper limit of these fluctuation amplitudes to be 
10~ 5 . The magnitude of the tensor perturbations depends upon the details of inflation potential [3(3] and is expected 
to be an order of magnitude smaller in amplitude than scalar perturbations. Therefore, we can set A x < 10 -6 . 
Hence, at T ~ 10 13 GeV with a very small amplitude of fluctuation such that A x ~ 10~ 19 the gravitational potential 
comes out to be Bq ~ 1CU 2 eV which is only an order of magnitude higher than the neutrino masses. Therefore, the 
corresponding oscillation probability and length are given by 



Vfi N = 0.999 shr 2 
AttE 



7 x 10~ 35 eV 2 sec 
AE 



A 



fIN 



7 x 10" 4 eV 2 



(99) 



which are an order of magnitude higher than that without gravitational effect. At the maximum possible A x , 
B ~ 10 11 eV >> m e , m„, m efi resulting the oscillation probability and length 



VfiN = sin 2 

XfIN 



1.4 x 10~ 22 eV 2 sec 



AE 



AttE 



1.4 x 10- 2 eVB ' 



(100) 



F. Implications 

The flavor mixing and hence oscillation is influenced by gravity. This may alter the relative abundance of different 
neutrino flavors in universe. Around a primordial black hole the maximum oscillation length of a thermal neutrino, 
from Eqs. (|88p and (f9T)]) . is obtained as 0.54 cm compared to 4.6 cm obtained from Eq. (|85p without gravity. 
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The oscillation probability increases as well significantly in early universe when we consider gravitational effect, as 
understood from section IV. D. At the neutrino decoupling era, while the probability increases for thermal neutrinos 
only 1.5 times, for TeV neutrinos it increases about two order of magnitude. It is also apparent from Eqs. (|84|). 

and (|94[) that at GUT era oscillation takes place vigorously. Therefore, production of muon neutrinos in early 
universe due to oscillation is expected to be much higher than that estimated without gravitational effect. 

At the inflationary era with primordial fluctuations, the oscillation length of, e.g., GeV neutrinos could vary from 
10 8 cm to 10 -5 cm depending on Bq, while size of universe is 10 3 cm. Therefore, the oscillation would be feasible. 
However, lifetime of the era itself is very small which may hinder significant oscillation. 

V. NEUTRINOLESS DOUBLE BETA DECAY 

It is generally believed that if neutrino is its own antiparticle, i.e. neutrino is a Majorana particle, then 
ncutrinolcss double beta decay may be observed. As CPT violating nature alters the mixing angle and masses in 
neutrino and antineutrino sector, it is interesting to see the effect of gravity on neutrinoless double beta decay rate. If 
one considers only one flavor, then gravity does not alter the decay rate. It was already shown [7| that CPT violating 
term has no effect on decay rate with only one flavor. Here we investigate the case for two flavor mixing. 

If we consider neutrino flavor mixing, without neutrino-antincutrino mixing, then one can express neutrino of 
different flavors as superposition of different mass eigenstates like (considering only electron and muon neutrino) 




(101) 



where U = J-(B° = 0), and fx and ]i are mass eigenstates with masses l\ and I2 respectively. In this case, the 
amplitude for neutrinoless double beta decay involving only electrons is |31j 

A cx J2kU^, i = 1,2. (102) 

i 

In presence of gravity the state \I/ is replaced by v as discussed in previous sections, and so the relation given by the 
Eq. (fTUTj) is modified to 




(103) 



Hence, the neutrinoless double beta decay amplitude is given by 

-4(1,2) oc M (li3) cos 2 0(1,2) + Af (2i4 )Sin 2 (/) (l! 2). (104) 
From the expressions of M-s and </>-s from Eqs. (|5Tj) and (|52[) and using Eqs. (|4"5|) and (j4T)|) we obtain 

Aoc ^Bl+ml. (105) 

Hence, under the effect of gravity the neutrinoless double beta decay amplitude is modified as the mass of the flavor 
state is modified. In absence of gravity the amplitude is simply proportional to m e . 

The double beta decay is involved with weak interactions, conventionally it has nothing to do gravity. Therefore, 
when we compute the decay rate with ^P-s, which are the states produced by pure charge-current interaction, it 
remains unchanged in presence of gravity. However, if we assume that under strong gravity ^-s modify to v-s and 
recalculate the rate, then it modifies by means of Bq. 

VI. SUMMARY AND DISCUSSION 

The idea of neutrino flavor oscillation is well established both theoretically and experimentally. However, there are 
some anomalies in different types of flavor oscillation data. For example, the LSND data provides mixing angle and 
neutrino masses which are inconsistent with solar and atmospheric neutrino oscillation data. This anomaly can be 
removed by introducing fourth family of sterile neutrino or CPT violation in neutrino sector. In a very past work, 
neutrino-antineutrino oscillation was discussed in analogy with K° — K° oscillation with maximal mixing [32j . In the 
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present paper, we have shown that in presence of gravity neutrino-antineutrino oscillation does occur. Neutrino and 
antincutrino interact differently with gravity through CPT violation. Under this condition if the Lagrangian has a 
lepton number violating Majorana mass, then neutrino-antineutrino mixing and oscillation occur. Without gravity, 
Majorana mass mixes neutrino and antineutrino states in equal ratio. Even if there is no Majorana mass, gravity 
splits neutrino and antineutrino in two different mass eigenstates. The gravitational scalar potential -Bo modifies the 
mass of the neutrino states. This Bq behaves like a lepton number conserving mass. The mass eigenstates, produced 
as a result of mixing under gravity, act as modified neutrino states. With these states we have studied two flavor 
mixing and corresponding oscillation generating four mass eigenstates. 

It is quite possible that the neutrino number what we see today is carried out off the neutrino decoupling era and the 
oscillation probability, strongly influenced by gravity, at that era should determine the relative abundance of today's 
flavors. We, in fact, have shown that the probability of conversion of electron neutrino to muon neutrino has greatly 
been enhanced, upto a few order of magnitude, with inclusion of gravity effect. On the other hand, even in present 
age the relative abundance of muon neutrino may be increased around primordial black holes where gravitational field 
is not negligible. 

It has been suggested that during core-collapse of a massive star, part of the infalling material goes into orbit 
around the compact core to form a hot, dense, centrifugally supported accretion disk whose evolution is strongly 
influenced by neutrino interactions. Under a wide range of conditions, this neutrino-dominated accretion flow will 
help to produce a successful supernova explosion (see, e.g. 34]) It will be interesting to study the effect of gravity on 
oscillation in such systems and then on related supernova explosion. 

In early, it was shown by pure quantum field theoretical consideration of neutrino flavor mixing that in infinite 
volume limit, the vacuum expectation value is not invariant under the transformation of flavor mixing (e.g. [35l. l36l| ) . 
This is the consequence of unitary inequivalance of the flavor and mass vacua. It has been shown that the vacuum 
structure with neutrino mixing has non-zero contribution to vacuum energy [37j | . It is well known that one of the 
interpretations of cosmological constant is linked with the density (p) and pressure (p) of vacuum with equation of 
state p = —p. The non-zero value of cosmological constant and hence the presence of such a vacuum energy density, 
which is often termed as dark energy, is needed to explain the observed acceleration of present Universe. Thus, the 
non-zero flavor vacuum energy density can be interpreted as a contribution to cosmological constant and hence to dark 
energy (35| . It has been shown 35 1 that the vacuum energy, the flavor vacuum expectation value of energy momentum 
tensor element Too, is non-zero and proportional to square of one of the Bogoliubov coefficients Vk associated with the 
flavor creation and annihilation operator. In ultrarelativistic limit |Vfc| 2 is proportional to (mi — TO2) 2 , the square of 
the difference of mass eigenvalues. Moreover, this vacuum expectation value is proportional to sin 2 <f>, where <fi is the 
mixing angle. Hence this contribution depends on the specific nature of the mixing. It is to be noted that in the limit 
Vk. — > 0, which is the case of traditional phenomenological mixing, the vacuum energy vanishes. All these calculations 
are in flat space. If we consider any curved spacetime, then the vacuum expectation value changes accordingly. Hence 
the contribution to vacuum density and then to dark energy depends on specific nature of mixing and background 
metric. In the present work, we have shown that background curvature i.e the gravitational field affects the mixing. 
Hence, in this case, we expect that gravity will play in both the ways - by affecting the mixing and by giving a curved 
background, to affect the contribution of dark energy. Therefore, it will be good to see how the present results can 
account for observed proportion of dark energy under the pure quantum field theoretic consideration. 

With this new mixing scheme, neutrinoless double beta decay rate has been revisited. With one neutrino flavor, 
it was shown [7] earlier that CPT violating term has no effect on decay rate. In this spirit we have calculated the 
amplitude of neutrinoless double beta decay considering neutrino-antineutrino admixture with two neutrino flavor 
mixing and still have found that the decay amplitude is proportional to the mass of the flavor state. However, as the 
mass of the flavor state, which is considered to be modified due to gravity, is now different from that in absence of 
gravity, the decay amplitude also differs from that in absence of gravity. Hence, the CPT violating gravity has an 
effect on the decay rate. 
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